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PREFACE 


Tables 1, 2 and parts of 3 and 4 of this publication were com- 
puted on the IBM 704 computer at the Berkeley campus of the Univer - 
sity of California. Because of the value of these tables in the applica- 
tion of cnoidal wave theory, the Council on Wave Research has prepared 
this publication. Although only the equations based on the work of 
Korteweg and de Vries are tabulated, a Similar set of tables derived 
from the work of Keulegan and Patterson have also been computed. This 
latter set of tables is not available in published form at this time. It is 
cautioned that these tables are unchecked. 

The authors wish to thank P. C. Klingeman for his suggestions 
during the preparation of the computer program, and R. Tokerud and 
Barbara A. Hyde for preparation of the illustrations. 


Frank D.Masch 
R..L. Wiegel 


University of California 
Berkeley, California 
March 1961 
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INTRODUCTION 


In 1895, Korteweg and de Vries developed the cnoidal wave 


theory and showed that it describes a class of permanent type long waves 


of finite amplitude. This theory which involves the Jacobian elliptic func- 


tions and gives the solitary wave and Sinusoidal wave as its two Lignans 
cases, appears to describe the progression of periodic waves into water 
whose depth is less than about 1 /8 or 1/10 the wave length more accur- 
ately than does the well-known theory for Stoke's waves. As indicated 


by Keulegan (1950), the theory of cnoidal waves is based on the assump- 


tion that the square of the inclindation of the water surface is small in re- 


lation to unity whereas for Stoke's waves to be valid in shallow water , 
the wave height must be very small. This restriction for Stoke's waves 
would appear to make the application of this theory somewhat unrealistic 
for shallow water waves. 

The cnoidal wave theory has been studied more recently by 
Keulegan and Patterson (1949), Keller (1948), Benjamin and Lighthill 
(1945), Iwasa (1955), Littman (1957), Wehausen and Rey 960), and 
Laitone (1960; 1961). Despite these rather extensive studies, compara- 
tively little use has been. made of the theory. In 1960, Wiegel summar- 
ized much of the existing work on cnoidal waves and presented the lead- 
ing results ina more usable form. However solutions of Wiegel's equa- 
tion for such wave characteristics as wave height, wave length, period, 


and celerity involve the use of graphs of many of the functions. To 


simplify further the application of the theory, this publication presents 
in tabular form some of the functions common to cnoidal waves. 

Wave characteristics have been tabulated from the equations ob- 
tained by Wiegel and based on the work of Korteweg and de Vries. The 
tables have been limited to the approximate range of validity based on 
the work of Laitone (1961) as shown in Fig. 1.. Also included are ible 
of the Jacobian elliptic functions. These functions have been extended 
over the range necessary for cnoidal wave application. Further dis- 
cussion of the format and use of the tabulated functions is given prior to 


the tables. 


SUMMARY OF CNOIDAL WAVE EQUATIONS 


The symbols used in the equations given below are defined as fol- 


lows: 

C wave speed 
d still water depth 
E(k) complete elliptic integral of the second kind 
H wave height 
k modulus of elliptic integral 
K(k) complete elliptic integral of the first kind 

LL wave length 
a wave period 
u horizontal component of water particle velocity 
u incomplete elliptic integral of the first kind 
Vv vertical component of water particle velocity 
x horizontal coordinate 
y vertical coordinate measured from the bottom 
Ys vertical distance from bottom to wave surface 
Yt vertical distance from bottom to wave trough 


The following equations for the wave characteristics are those ob- 


tained by Wiegel and are based on the theory of Korteweg and de Vries. 
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In dimensionless form, the general expression for the wave length is 


WSR oie 2 
+: = [k K(K)] 


q3 GL) 


From Stoke's second definition of wave velocity which is the velocity of 
propagation of the wave form when the horizontal momentum is reduced 


to zero by the addition of uniform motion, the wave celerity is given by 


petal BH 
Z K(k) 


As cnoidal waves are periodic and of permanent type, .1e period is 


k K (k) 
H 


(Ua Oa er iy Bes 
Decree Beto dk2 E Ko] (3) 


Eqs. (1), (2), and (3) involve the complete elliptic integrals of the 


first and second kind. Although the values of these integrals have been 
tabulated over the range k2 = 1 - io-& by Kaplan (1946; 1948), Hayashi 
(1930; 1933), and Airey (1935), the application of cnoidal wave theory re- 
quires that the range be extended to k2 = 1-10°*9, This was-done by 


Wiegel using the following series from Jahnke and Emde (1945). 
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where k'= \V(1 -k2) and A is the natural logarithm of 4/k'. 


Cnoidal wave profiles are determined from the equation 
“ 2 CiaD egies 
Ys = yt + Hen ppKc L pT) ; (4) 


where cn ( ) is the Jacobian elliptic function, x/L is the same as u/2K(k) 
for the cn@() function, and U is the incomplete elliptic integral of the 
first kind. u has been used rather than the more commonly accepted 
symbol u to avoid confusion with the horizontal component of water 
particle velocity. Values of these functions which are often written as 
cn(t |x), have been tabulated over a limited range of k? by Milne- 
Thompson (1950), Schuler and Gabelein (1955), and Spenceley and 
Spenceley (1947). The range of the Jacobian elliptic functions has also 


been extended by Wiegel with the following equation from Milne Thompson 


4 tanh @ sech u (sinh u cosh JU - @) 


en (7 |) = sech u - 4 k! 
Typical wave profiles computed from Eq. (4) are shown in Fig. 2. 


The horizontal and vertical components of water particle velocity 


are determined from the following equations of Keulegan and Patterson 
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celerations are 
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The use of the above equations for velocities and accelerations requires. 
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following identities between the Jacobian elliptic functions: 


T|k = 2K(k) + (11) 
sn2(u | k) = 1 - cn#(aj k) (12) 
dn2(@ {k) 21 -k? [1 - on2(t|k) | (13) 


Examples of horizontal and vertical components of water particle 


velocities and local accelerations are shown in Figs. 3 - 6. 


DESCRIPTION OF TABLES 

It is the purpose of this section to give a brief discussion of the 
format and the use of each of the tables. Also included are graphs show- 
ing the relations between the various wave characteristics and some com- 
binations of elliptical functions. Since the use of the enue is probably 
best illustrated by an example, a simple problem is worked out at the 
end of this section. 

As previously mentioned, the application of cnoidal wave theory re- 
quires that the range of k2 be extended to k2 = 1 - 10749. Values of k@ 
up to k? = 1 - 107° are tabulated in decimal Rev at SE. 5. ML 107° is 
written 0.99999. For k2 greater than 1 - Lo"? , the exponential notation 
is used. 

As the increments used in computing the tables are small, a linear 
interpolation for intermediate values should be satisfactory. If however, 
additional accuracy is desired, a second order scheme of interpolation 


should be employed. 
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Trevouley Ibe 

This table lists values of L2H/d®, k, k2, and K(k). Values of 
L7H/d? were computed from Eq. (1) given in the summary of cnoidal wave 
equations. Table 1 is used by entering with either a known value of k2 as 
determined from a different table, or with L2H a3 from which the value of 
k? is obtained. The graphical relation between the parameters listed in 
“aole lace Shown an Eis. 7. 
Mable 2. 

Table 2.is a series of individual tables in which values of T\/g/d, 
C*/ ed, k, k*, K(k), E(k), and E(k)/K(k) are tabulated for multiples of the 
parameter, H/d. The wave characteristics, T\g/d and C2/g¢d, were com - 
puted from Eqs. (3) and (2) respectively. The tables have been limited in 


accordance with Fig. 1, and the range of H/d is from 0.01 to 0.78 in incre- 


ments of 0.01. These are the values useful in the valid range for the cnoidal 


theory. The remaining portions of the table are on file at the University of 
California's Hydraulic Engineering Laboratory. Plots of the wave charac- 
teristics are given in Fies.8, 9, and 10. -fable.2 is used by entering 

with H/d and either TYVg/d or C2/gd from which k2 is determined. The 
table may also be entered with H/d and a value of k# in which it is possible 
to read off the values of TVg/d and C2/ed. A third way in which Table 2 
can be used is to enter with k2 and either T g/d or C*/2d .) In this anenner, 
the value of H/d is obtained. 

Table 3. 


Table 3 gives the Jacobian elliptic functions for multiples of k2 and 
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7. These functions are computed for u over the range of 0.00 to 8.00 in 
increments of 0.05. Eqs.(11), (12), and (13) were used in the computa- 
tion of this table. Also included in Table 3 is the parameter u/2K(k) and 
the product of the Jacobian elliptic functions. 

Table-4. 

As there is very little change in the value of the Jacobian elliptic 
functions for 1 - 1076 < ke —< 1 - OP ee Table 5 assumes these functions 
to be constant over this range. (A calculation of the sn, cn, and dn func- 
tions at k? = 1 - 10-40 shows no change in the value of the functions through 
the five decimal places used in the table.) Therefore, Table 4 effectively 
gives only the values of u/2K(k) for multiples of u, and for k@ in the range 
of 1 - 10-6 to1 - 10°*9. This table is used by entering with the values of 
kK? and u from which the constant Jacobian elliptic functions, their product, 
and u/2K(k) are obtained. 

Tabless . 

In order to simplify the use of the tables of sn, cn, and dn avyail- 
able in Milne-Thomson and in Table 4, Table 5 has been prepared with 
u/2K (i.e., x/L) tabulated as a function of u and k?. 

Ulustrative Example 

The following example will illustrate the use of the tables for so- 
lution of the wave characteristics. 

Given: Wave height = 0.5 ft. 


Wave length = 25 ft. 
Water depth = 2.0 ft. 


Find: Wave period 
Wave velocity 


IE 


Solution: 
L2H/d = 39.06 
iil = O28 
From Table 1, k? = .936 ; 
From Table 2, f 
a Visicl = Po ; 

2 

Gy fede 1,06 


Therefore, | i 


ae 


3.02 sec. 


C 


1 
8.28 ft. per sec. 
4 
q 


ILLUSTRATIONS OF SOME USEFUL COMBINATIONS 
OF FUNCTIONS 


Some combinations of sn, cn, dn functions together with the 


modulus kK? occur in equations for water particle velocities and accel- 


erations; some of these combinations have been plotted in Figs. 11-17. 
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